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g i Abstract 

The high-temperature plasma above the electroweak scale ~ 100 GeV may have contained 
a primordial hypercharge magnetic field whose anomalous coupling to the fermions induces 
a transformation of the hypermagnetic energy density into fermionic number. In order to 
describe this process, we generalize the ordinary magnetohydrodynamical equations to the 
anomalous case. We show that a not completely homogeneous hypermagnetic background 
induces fermion number fluctuations, which can be expressed in terms of a generic hyper- 

X ■ 

magnetic field configuration. We argue that, depending upon the various particle physics 
parameters involved in our estimate (electron Yukawa coupling, strength of the electroweak 
phase transition) and upon the hypermagnetic energy spectrum, sizeable matter-antimatter 
fluctuations can be generated in the plasma. These fluctuations may modify the predictions 
of the standard Big Bang nucleosynthesis (BBN). We derive constraints on the magnetic 
fields from the requirement that the homogeneous BBN is not changed. We analyse the 
influence of primordial magnetic fields on the electroweak phase transition and show that 
some specific configurations of the magnetic field may be converted into net baryon number 
at the electroweak scale. 
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1 Introduction 



There are no compelling reasons why magnetic fields should not have been present in the 
early Universe. Moreover, it can be argued that the existence of some magnetic fields at high 
temperatures is a quite natural phenomenon. Indeed, the presence of large scales magnetic 
fields in our observed Universe is a well established experimental fact. Since their first 
evidence in diffuse astrophysical plasmas beyond the solar corona fl], ||, magnetic fields have 
been detected in our galaxy and in our local group through Zeeman splitting and through 
Faraday rotation measurements of linearly polarized radio waves. The Milky Way possesses 
a magnetic field whose strength is of the order of the microgauss corresponding to an energy 
density roughly comparable with the energy density today stored in the Cosmic Microwave 
Background Radiation (CMBR) energy spectrum peaked around a frequency of 30 GHz. 
Faraday rotation measurements of radio waves from extra-galactic sources also suggest that 
various spiral galaxies are endowed with magnetic fields whose intensities are of the same 
order as that of the Milky Way ||. The existence of magnetic fields at even larger scales 
(intergalactic scale, present horizon scale, etc.) cannot be excluded, but it is still quite 
debatable since, in principle, dispersion measurements (which estimate the electron density 
along the line of sight) cannot be applied in the intergalactic medium because of the absence 
of pulsar signals [|J. 

If the existing galactic magnetic field is naively blue-shifted to earlier epochs, one then 
finds that the Universe was always filled by a magnetic field of a considerable amplitude, with 
energy density of the order of the energy density of 7 quanta. Of course, this consideration 
does not take into account different physical processes operating at the galactic scale, such 
as the dynamo mechanism [0, |J] or the anisotropic collapse mechanism ||, which change 
considerably the naive scaling law \H\ ~ T 2 . In any case, it is widely believed that some 
seed fields of primordial origin are necessary for the successful generation of the galactic 
magnetic fields 0. 

Looking at this problem from a more theoretical side, there are several mechanisms that 
may successfully generate large enough magnetic seeds coherent on different scales. Mag- 
netic seeds can be produced either during a first-order quark-hadron phase transition j/], § 
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or during the electroweak phase transition || |10|, |TTj, |I2| . Recently it was also suggested 
that a primordial asymmetry encoded in the right electron number can be converted in a 
quite large hypercharge seed during the symmetric phase of the electroweak theory 
The seeds could also be the result of the parametric amplification of the quantum mechan- 
ical (vacuum) fluctuations of some primordial gauge field, in the same way as in General 
Relativity the quantum mechanical fluctuations of the tensor modes of the metric can be 
amplified, producing, ultimately, a stochastic gravity- waves background JTJ]]. The essential 
ingredient of the large-scale magnetic field generation is the breaking of conformal invariance 



in the coupling of the electromagnetic field to gravity JToJ, jig] . Reasonable seeds could also 
be produced if the inflaton is coupled to the Maxwell term in a chaotic inflationary scenario 
In the string theory low-energy effective action, the dilaton field provides a unified 
value of the gravitational and gauge coupling at the string scale and it naturally breaks the 
scale invariance of the electromagnetic and gauge couplings (also in four dimensions) with- 
out providing a gravitational mass for the photon. Sizeable seeds, coherent over the galactic 



scale, can be generated |18], |19 



The possible existence of magnetic fields in the early Universe has a number of interesting 
cosmological implications. For example, magnetic fields at small scales may influence the 
Big Bang nucleosynthesis (BBN) and change the primordial abundances of the light elements 
|20| by changing the rate of the Universe expansion at the corresponding time. The success 
of the standard BBN scenario can provide an interesting set of bounds on the intensity of 
the magnetic fields at that epoch . 

Long-range stochastic magnetic fields that possibly existed at the decoupling epoch might 
have induced anisotropies in the microwave sky [T^J. The existence of a completely homoge- 
neous field coherent over the horizon at the present epoch can be interestingly constrained by 
the COBE observations |21], p2]] . Conversely, if the CMBR is linearly polarized, its polariza- 
tion plane could have been rotated by the presence of a sufficiently energetic magnetic field 
coherent over the horizon size at the decoupling epoch |23| . Faraday rotation measurements 
applied to the galactic (synchrotron) emission can also provide interesting constraints |24 
on large-scale magnetic fields (even though these are coherent over scales smaller than the 
present horizon). 
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In this paper we address the question of whether there can be any cosmological conse- 
quences from the fact that magnetic fields existed prior to the electroweak phase transition, 
when the background temperature was T > T c ~ 100 GeV. At these temperatures the elec- 
troweak symmetry is restored and the magnetic field is replaced by the hypermagnetic one. 
The hypercharge field, unlike the ordinary magnetic field, has an anomalous coupling to the 
fermions. This fact will play a crucial role in our considerations. The origin of primordial 
hypermagnetic fields is not essential for us and consequently we simply assume that they 
were generated by some mechanism before the electroweak (EW) phase transition. 

We will show that, depending on the particle physics model and on the initial spectrum 
of the primordial magnetic fields, quite large fluctuations of the baryon and lepton numbers 
may be generated. These fluctuations can survive until the onset of BBN and create unusual 
initial conditions for the calculation of the light element abundances. In particular a natural 
outcome of our considerations is a non-uniform distribution of baryon number, not necessar- 
ily positive-definite. Matter-antimatter domains are then possible. The requirement that 
these fluctuations are small enough in order not to conflict with the predictions of the stan- 
dard homogeneous BBN allows us to put quite a strong constraints on the spectrum of the 
primordial magnetic fields. Moreover we will argue that the primordial magnetic field may 
change the nature of the electroweak phase transition. Finally, the existence of primordial 
fields with some specific topological structure may result in the production of the net baryon 
number at the electroweak scale. 

The plan of the paper is the following. In Section 2 we derive our basic equations. 
For an ordinary electromagnetic plasma, it is fairly well established that the evolution of 
the magnetic fields can be described using the magnetohydrodynamical (MHD) equations 
H, f|, |25fl . In the case of a high-temperature electroweak plasma the MHD equations have 
to be generalized by taking into account anomaly effects. In the new AMHD equations 
(anomalous MHD equations) the magnetic hypercharge fields turn out to be coupled to 
the fermionic number density. As a consequence, the evolution equations of the anomalous 
charge densities acquire a magnetic source term. In Section 3 we describe an approximate 
solution of AMHD equations. We show that anomalous coupling gives rise to instabilities, 
allowing the conversion of the energy sitting in the fermionic degrees of freedom into magnetic 
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hypercharge fields and vice versa. This phenomenon is completely absent in ordinary MHD. 
The presence of the fermionic number density produces a kind of "ohmic" current, which is 
parallel to the magnetic hypercharge field. Also this phenomenon is quite new if compared 
with the ordinary MHD equations (though something vaguely similar can happen in the 
context of the dynamo mechanism || in a parity- non- invariant turbulent fluid). In Section 
4 we will apply our results to the case of stochastic hypercharge field backgrounds, whereas 
in Section 5 we focus our attention on the possible phenomenological implications of our 
considerations for BBN. Different bounds on primordial magnetic fields will also be analysed 
in the light of our results. We will also discuss the dependence of the magnitude of the 
baryon number fluctuations upon various particle physics parameters, which can appear in 
the extensions of the minimal standard model (MSM). In Section 6 we discuss the influence 
of the primordial hypermagnetic fields on the dynamics of the electroweak phase transition, 
and we will show how some specific configurations of the hypermagnetic fields may create a 
net non- zero amount of baryons. Section 7 contains our concluding remarks. 

Part of the results of this paper has been already presented (in a more compact form) in 



(see also the closely related paper ]T3| where transformation of a finite number density 
of right electrons into hypermagnetic fields has been considered). 



2 Basic equations 
2.1 Preliminaries 

Let us start from some qualitative considerations. A unique property of "unbroken" U(l) 
gauge interaction is the absence of mass of its corresponding vector particle. Static "mag- 
netic" fields are never screened (in the absence of monopoles) and thus homogeneous fields 
can survive in the plasma for infinitely long times. On the contrary, electric fields quickly 
decay because of the finite plasma conductivity a within a time scale ~ 1/cr. Then long- 
ranged non-Abelian magnetic fields (corresponding to, e.g. the colour SU(3) or weak SU(2) 
groups) cannot exist because at high temperatures the non-Abelian interactions induce a 
"magnetic" mass gap ~ g 2 T. Also the non-Abelian electric fields decay because of the finite 
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value of the conductivity as it occurs for Abelian electric fields. Therefore, the only long 
scale field that can exist in the plasma for enough time must be associated with some Abelian 
U(l) group. This statement, valid to all orders in perturbation theory, has been confirmed 
non-perturbatively for the electroweak theory by recent lattice studies in p7 |. Under normal 
conditions (i.e. small temperatures and small densities of the different fermionic charges) 
the SU(2)xU(l)y symmetry is "broken" down to U(1)em, the massless field corresponding 
to U(l) em is the ordinary photon and the only long-lived field in the plasma is the ordinary 
magnetic one. At sufficiently high temperatures, T > T c , the SU(2)xU(l)y symmetry is 
"restored", and non-screened vector modes correspond to the U(l)y hypercharge group. 
Hence, if primordial fields existed at T > T c , they did correspond to hypercharge rather than 
to U(1)em- 

There are essential differences between the interactions of magnetic fields and the ones 
of hypermagnetic fields with matter. The ordinary electromagnetic field has a vector-like 
coupling to the fermions, while the coupling of the hypercharge fields is chiral. Thus, if 
hyperelectric (£y) and hypermagnetic (7iy) fields are present simultaneously, they cause a 
variation of the fermionic number according to the anomaly equation, d^j^ ~ f^T^y ■ £y 
(here g' the hypercharge gauge coupling constant). Now, the presence of non-homogeneous 
hypermagnetic fields in the EW plasma with (hyper) conductivity a c always implies the 
existence of a related electric field, Ey ~ ^-V x Hy Since for a general stochastic magnetic 
background {{H Y ■ V x H Y ) 2 ) ^ 0, the non-uniform hypermagnetic field may absorb or 
release fermions and produce, ultimately, baryon and lepton density perturbations because 
of the anomaly equation. 

The behaviour of cold fermionic matter with non-zero anomalous Abelian charges was 



considered in where it was pointed out that the anomalous fermionic matter is unstable 
against the creation of Abelian gauge field with non-zero Chern-Simons number, which eats 



up fermions because of the anomaly. It was suggested in |U| that the right electron number 
density may be converted to the hypercharge field because of a similar effect. The main 
aim of this paper is the study of the opposite situation, namely we want to understand how 
hypercharge fields may be converted into fermions in a hot environment. 
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2.2 MHD equations for ordinary plasmas 

During the symmetric phase of the electroweak theory the evolution of the background 
geometry is dominated by radiation. The first assumption we will make is that prior to 
T c ~ 100 GeV the geometry may be described by a conformally flat metric of the FRW type, 
whose line element is: 

ds 2 = g af3 dx a dx p = a(r) 2 (dr 2 - dx 2 ), a(r) ~ r (2.1) 

(r is the conformal time coordinate related to the cosmic time coordinate t as dt = a(r)dr). 
We will also assume that the radiation-dominated stage started much before the electroweak 
epoch at temperatures T > T c . 

The Weyl invariance of the ordinary Maxwell equations in a conformally flat FRW back- 
ground geometry implies that the MHD equations in the metric (|2.1| ) can be written |29) 
as: 

dH ^ - BE -> 

— = -VxE, +j = VxH, 

or or 

V-H = 0, V-E = 0, 

V-J = 0, J = a(E + vxH) (2.2) 

(H = a 2 H, E = a 2 £; A = aA; J = a 3 ]; a = a c a; H = V x A; £ = —^A; H, £, A, 
j, <r c are the flat-space quantities whereas H, E, A, J, a are the curved-space ones; v is 
the bulk velocity of the plasma). We assumed that the plasma is locally electrically neutral 
(V • E = 0) over length scales larger than the Debye radius. We notice that the spatial 
gradients used in Eq. ( |2.2|) are defined according to the metric ( |2.1| ). 

There are several approximations in which the above equations can be studied. One is 
the so-called ideal (or superconducting) approximation and the other is the real (or resistive) 
case 



25, 30 



In the ideal case a" 1 = and, from Ohm's law we can immediately deduce that the 
electric field is orthogonal to the magnetic one and it is also orthogonal to the bulk velocity 
of the plasma: 

E~-vxH. (2.3) 
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Two important theorems of the ideal MHD, which follow from Eq. ( |2.3|) , are the conservation 
of the magnetic flux: 

$= f H-dE (2.4) 



and of the magnetic helicity (Chern-Simons number) || p0[ : 

Ncs = —I d 3 x H ■ A, a em = (2.5) 

where <i£ is a closed surface in the plasma; the volume integral is performed over a magnetic 
flux tube. 

If, on the contrary, the effect of the finite value of the conductivity is taken into account 
(a -1 <C 1) and the resistive Ohm law is employed, then the induced electric field is not 
exactly orthogonal to the magnetic one and the conservation laws of the ideal MHD are 
corrected (in the resistive approximation) by an expansion in powers of the resistivity, which 
can be explicitly computed and which we report at the lowest order in <j~ l : 

d _ 1 



— $ = — / V x (V x H) ■ tffi, 
dr a Jt, 



d . Ncs = -^l I d'xH-VxH. (2.6) 



dr ncr jv 

According to Eqs. ( |2.6j ) the magnetic flux lines evolve glued to the hypercharged plasma 
element; also the sum of the link and twist number of the magnetic flux tubes is always the 
same all along the time evolution, only provided that a -1 = 0. 

The same dynamical information encoded in the magnetic flux conservation is also con- 
tained in the magnetic diffusivity equation, which can be derived according to Eq. (|2[ 



an i 

^L = Vx(vxH) + -V 2 H. (2.7) 
or o 

The ratio of the two terms on the r.h.s. defines the magnetic Reynolds number 

a\VxvxH\ 

R ~ — ~ v I a. 2.8) 

\V 2 H\ 

If R <C 1 (for a given length scale I) the flux lines of the magnetic field will diffuse through 
the plasma. If R 3> 1 the flux lines of the magnetic field will be frozen into the plasma 
element. From the magnetic diffusivity equation (|2.6|) it is possible to derive the typical 



structure of the dynamo term by carefully averaging over the velocity field according to the 
procedure outlined in [Q, |3J]. By assuming that the motion of the fluid is random and has 
zero mean velocity, it is possible to average over the ensemble of the possible velocity fields. 
In more physical terms this averaging procedure of Eq. (|2.7|) is equivalent to averaging over 
scales and times exceeding the characteristic correlation scale and time tq of the velocity 
field. This procedure assumes that the correlation scale of the magnetic field is much larger 
than the correlation scale of the velocity field. In this approximation the magnetic diffusivity 
equation can be written as: 

^ = a(Vxff) + -V 2 H (2.9) 

OT O 

(a = — 2 j(v-'V xv) is the so-called dynamo term, which vanishes in the absence of vorticity; in 
this equation H is the magnetic field averaged over times larger than tq, which is the typical 
correlation time of the velocity field). We can clearly see that the crucial requirement for the 
all averaging procedure we described is that the turbulent velocity field has to be "globally" 
non- mirror- symmetric. It is interesting to point out that the dynamo term in Eq. (|2.9| ) 
has a simple electrodynamical meaning, namely, it can be interpreted as a mean ohmic 
current directed along the magnetic field: 

J = -aH. (2.10) 

This equation tells us that an ensemble of screw-like vortices with zero mean helicity is able 
to generate loops in the magnetic flux tubes in a plane orthogonal to the one of the original 
field. This observation will be of some related interest for the physical interpretation of the 
results we are going to present in the following paragraph. We finally notice that if the 
velocity field is parity-invariant (i.e. no vorticity for scales comparable with the correlation 
length of the magnetic field), then the dynamics of the infrared modes is decoupled from the 
velocity field since, over those scales, a = 0. 

2.3 AMHD equations for electroweak plasmas 

The electroweak plasma in complete thermal equilibrium at a temperature T can be char- 
acterized by rif chemical potentials /ij, i — 1, ...,rif corresponding to the exactly conserved 



S 



global charges 



Ni = L { 



(2.11) 



(Li is the lepton number of the z-th generation, B is the baryon number, and n/ is the number 
of fermionic generations). One should also introduce a chemical potential fiy corresponding 
to weak hypercharge, but its value is fixed from the requirement of the hypercharge neutrality 
of the plasma, (Y) = 0. 

We want to study this plasma slightly out of thermal equilibrium; in particular, we want 
to see what happens with a non-uniform distribution of the hypermagnetic field. Because of 
the anomaly, and thanks to the arguments illustrated in Section 2.1, this field is coupled to 
the fermionic number densities. In principle, different chemical potentials can be assigned to 
all the fermionic degrees of freedom of the electroweak theory (45 if n/ = 3) and the coupled 
system of Boltzmann-type equations for these chemical potentials and the hypercharge fields 
may be written. Since we are interested in the slow processes in the plasma, this is not 
necessary. If the coupling, corresponding to some slow process, is switched off, then the 
electroweak theory acquires an extra conserved charge and a further chemical potential 
should be added to the system given in Eq. fl2.11|) . 

An interesting observation (which turns out to be quite important in our context) has 
been made in |HJ], where it was noticed that perturbative reactions with right-handed elec- 
tron chirality flip are out of thermal equilibrium at temperatures higher than some temper- 
ature Tr.| Thus, the number of right electrons is perturbatively conserved at temperatures 
T > Tr and the chemical potential hr can be introduced for it. On the other hand, this 
charge is not conserved because of the Abelian anomaly, 



and it is therefore coupled to the hypermagnetic field. Here y and y are, respectively, the 
Uy(l) hypercharge field strengths and their duals, g' is the associated gauge coupling and 
yji = —2 is the hypercharge of the right electron. 

2 This temperature depends on the particle physics model, see also the discussion reported in Section 5. 
In the MSM T R ~ 80 TeV |l|. 
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y^\ 



(2.12) 
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Now we are ready to derive the anomalous MHD equations |T3|, |26|]. The effective La- 



grangian density describing the dynamics of the gauge fields at finite fermionic density is 



[32 
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C Y ,e R = -\^Y aP Y al3 - ^g~J a Y a + fxe ijk Y^Y k , fi = J-^ R (2.13) 

(g is the determinant of the metric defined in ( |2.1| ); Y af3 = Vr a Yg]; V Q is the covariant 
derivative with respect to the metric fl2.1|) [notice that in the metric ( |2.1|) V[ a Y/3] = <9[a^/3]]; g' 
is the Abelian coupling constant). The first and the last terms in Eq. ( |2.13|) are nothing but 



the curved space generalization of the flat-space effective Lagrangian for the hypercharge 
fields at finite fermion density |32j, J a is the ohmic current. The equations of motion for the 



hyperelectric and hypermagnetic fields are then 

dH OE n' 2 

—L = -V x E Y , "~ 7r~ + J Y = -^maHy + VxH Y , 
or or 7T 

V • Hy = 0, V ■ Ey = 0, 

V • Jy = 0, J Y = a[E Y + v x H Y ), o = a c a(r) (2.14) 

(E Y = a 2 £ Y , Hy = a 2 H Y , Hy = V xy,£ Y = -^S). For the EW theory the conductivity 
of plasma is a ~ 70T ||33|| . To the equations of motion of the hypercharge field ( |2.14j ) we 



have to add the evolution equation of the right electron chemical potential, which accounts 
for the anomalous and perturbative non-conservation of the right electron number density 
(n R ) : ^ 

^jf = -£-2?r ■ ~ r K - (2.15) 
where T is the perturbative chirality-changing rate, T = Tj^, n e ^ is the equilibrium value of 
the right electron number density, and the term proportional to Ey ■ Hy is the right electron 
anomaly contribution. 

Finally, the relationship between the right electron number density and the chemical 
potential must be specified. This relation depends upon the particle content of the theory, 
e.g. upon the number of fermionic generations, the number of Higgs doublets, etc. We 



will write it for the case of the MSM [34]: the specific coefficients change only slightly for 
other theories and do not have a significant impact on the results. For generality, we assume 
that the Universe is asymmetric not only with respect to the number of right electrons, but 
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also with respect to corresponding conserved charges defined in (|2.11 ) and compute all the 
relevant chemical potentials: 

2 2Ar .783 . 201 . 15,. 

2 2Ar . 201 1227 . 3 ^ lrT1 

Mi = —irN eff \ 6 R -\ 5\ H (5 2 + 5 3 )}T, 

45 //L 88 440 110 V ;J ' 

^ 2 45 //L 22 H 110 55 55 3J ' 

u 3 = — irN eff \ — 5 R -\ 5 1 + -\ 5 2 -\ 8 3 ]T. 

^ 45 //[ 22 R 110 55 55 3J ' 



9 97 11 

a Y = ±^N eff [—S R + —5 1 + —(5 2 + 5 3 )}T, (2.16) 
45 //L 88 440 110 V 1 



where 5i is the asymmetry of z-th conserved charge, 5r = n R /s (s is the entropy density) is 
the right electron asymmetry, N e ff = 106.75 is the effective number of relativistic degrees of 
freedom in the symmetric phase of the MSM. 

With the use of relations ( 2.16Q , Eq. ( [2.15| ) can be rewritten completely in terms of the 
right electron chemical potential: 

a dr 4vr 2 88 a 3 T 3 K ^ H J K J 

At finite hyperconductivity (in what we would call, in a MHD context, "resistive" ap- 



proximation) we have that from Eq. ( 2.14 ) the induced hyperelectric field is not exactly 
orthogonal to the hypermagnetic one and, moreover, an extra "fermionic" current comes in 
the game thanks to the fact that we are working at finite chemical potential. Therefore in 
our context the resistive Ohm law can be written as 

Ey = — -VX Ey ~ - ( — hrclHy + V X Ry ) -VX Hy, a' = %- . (2.18) 
a a \ ix J Arc 

In the bracket appearing in Eq. (|2.18| ) we can identify two different contributions. One is 
associated with the curl of the magnetic field. We will call this the MHD contribution, since it 
appears in the same way in ordinary plasmas. The other contribution contains the chemical 
potential and it is directly proportional to the magnetic field and to the chemical potential. 
This is a typical finite density effect. In fact the extra Ohmic current simply describes the 
possibility that the energy sitting in real fermionic degrees of freedom can be transferred 
to the hypermagnetic field. It may be of some interest to notice the analogy between the 
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first term of Eq. (|2.18 ) and the typical form of the ohmic current discussed in Eq. (|2.10 ) 



appearing in the context of the dynamo mechanism. In the latter case the presence of a 
current (proportional to the vorticity through the a dynamo term) was indicating that large 
length scales magnetic fields could grow by eating up fluid vortices. By inserting E Y obtained 
from the generalized Ohm law (|2.18 ) in the evolution equations fl2.14|) of the hypercharge 



fields, we obtain the generalized form of the magnetic diffusivity equation 

dHy = _ x f g \ + v x (i; x #) + -V 2 H Y . (2.19) 

or ixo v ' a 

In order to be consistent with our resistive approach we neglected terms containing time 

derivatives of the electric field, which are sub-leading provided the conductivity is finite. In 

our considerations we will also make a further simplification, namely we will assume that the 

EW plasma is (globally) parity- invariant and that, therefore, no global vorticity is present. 

Therefore, since the length scale of variation of the bulk velocity field is much shorter than 

the correlation distance of the hypermagnetic field, the infrared modes of the hypercharge 

will be practically unaffected by the velocity of the plasma, which will be neglected when 

the large-scale part of the hypercharge is concerned. This corresponds to the usual MHD 

treatment of a mirror symmetric plasma (see, e.g. Eq. (|2.9| ), when a = 0). 

Eqs. ( |2.19| ) and ( |2.17| ) form a set of AMHD equations for the hypercharge magnetic field 
and right electron chemical potential in the expanding Universe. 

An important property of Eqs. (|2.14|) and ( |2.19|) is that they are not conformally invari- 
ant. The conformal invariance of the ordinary Maxwell equations implies that the equations 
for the rescaled fields in curved space keep the same form also in flat space in terms of the 
non-rescaled fields provided the conformal time coordinate is employed in curved space and 
the cosmic time coordinate is employed in flat space. We can easily see that this is not the 
case of Eq. (|2.14j ) by writing our evolution equations in flat space: 

-rrf = -VxS Y , —f + j Y = ^ flR H Y + VxH Y , 

Ot Ot 7T 2 

V • Hy = o, V • By = 0, 

V-jy = 0, J Y = a c £ Y . (2.20) 

The lack of conformal invariance comes from the presence of the scale factor in front of the 



12 



right electron chemical potential in the evolution equation Q2.14|) for Hy- Clearly, the explicit 



breaking of conformal invariance is also reflected in the Ohm law and in the hypermagnetic 
diffusivity equation which, passing from curved to flat space, become 

£y = J —~-(-vrH y + V><Hy), (2.21) 

^ = -— V x U R H Y ) + lv 2 7? y . (2.22) 
Ol no v ' cr c 

In flat space the kinetic equation of the right electron chemical potential becomes instead: 

dn R _ g' 2 783 -> -* 

~df ~ -I^T^ £y ' Hy ~ TfXR - (2 ' 23) 

It is interesting to notice that the term containing the chemical potential in Eq. (|2.19| ) 
plays a role similar to that of the dynamo term, since it also produces an instability [ |D| . 
Its physical interpretation is actually quite simple. We could define, as in the case of the 
ordinary MHD, a generalized Reynolds number that measures the relative weight of the 
two terms on the r.h.s. of Eq. ([2.19|) . If the diffusion term dominates, then the flux of 
magnetic hypercharge will diffuse through the plasma. If on the contrary we are in the inertial 
range (where the diffusivity term is negligible) there are two possibilities. If the chemical 
potential is exactly zero then the hypermagnetic field will be frozen into the plasma element 
as required by magnetic flux conservation. However if we are at finite fermion density the 
energy density sitting in fermionic degrees of freedom may be transformed in infrared modes 
of the hypermagnetic field. 

3 Fermions from the hypercharge field 

3.1 T >T C 

In this section we are going to compute the relationship between hypermagnetic fields and 
induced fermionic chemical potential, at temperatures larger than the critical temperature 
of the electroweak phase transition T c . 

There is an important consequence of the resistive approximation. By using the Ohm 



law given by Eq. (2.21|) we can eliminate the hyperelectric field from the kinetic equation of 
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the right electrons and obtain 

9 ffiR\ g' 2 783 - -* . .fi R 



*\TJ OT88- W - VxHl '- (r + r «»T' (2 - 24) 
where 

_ 783 a' 2 \H Y \ 2 
Th ~ -22^ — - (2 ' 25) 
We notice immediately that the source term appearing in the r.h.s. of Eq. ( p.24| ) (and 
coming from the anomaly) is indeed strongly reminiscent of what we would call (in a MHD 
context) magnetic helicity. From Eq. (|2.24j ) one can see that the right-electron number is not 
conserved (even if T = 0) because of the Abelian anomaly, provided a non-zero hypermagnetic 
field is present (cf. Ref. P5| ). Eq. Q2.24| ) can be solved in the adiabatic approximation at 
temperatures T < Tr, when perturbative right electron chirality flip reactions are in thermal 
equilibrium. Neglecting the time derivative of the chemical potential, we get: 

fjR a' 783 H Y -VxH Y 

T ~ na c T 3 88 T + T n [ ' ' 

The solution (|2.26| ) can be inserted into Eq. ( |2.22|) for the magnetic field, which will become 
a partial (non-linear) differential equation containing only the hypermagnetic field. Thus, an 
inhomogeneous hypermagnetic field can produce a spatial variation in the chemical content 
of the plasma. In fact, according to Eq. (|2.16|) , spatial fluctuations in the right electron 
chemical potential will determine fluctuations not only in the right electron number density 
but also in the number densities associated with the other fermion asymmetries. Fluctuations 
in the number density of some species are frequently called isocurvature perturbations. There 
are actually two regimes where Eq. ( 2.26| ) can be analysed. The first one is the regime 
where V > T^. In this case the rate of right electron dilution is essentially determined by 
the perturbative processes, which can flip the chirality of the right electrons. In the opposite 
case (T^ > T) the rate is mainly due to the presence of the Abelian anomaly, and 

88tt 2 {Hy-VxH y \ n ^ . 



783#' 2 v \h y \ 2 J vr 



H 



In Eq. (|2.27|) the chirality-changing rate only comes in the correction. Moreover, since 
the hypermagnetic field intensity appears with the same power in the numerator and in 
the denominator, the right electron fluctuations are independent of the magnitude of the 



14 



hypermagnetic field fluctuations and fixed by their spatial distribution. It is interesting to 
notice that in Eq. (2.27) the actual value of the conductivity completely cancels and only 
appears in the correction. 



3.2 T <T C 

Now we are going to discuss what happens after the electroweak phase transition. Below T c 
the hypercharge magnetic fields are converted into ordinary magnetic fields. The latter do 
not have anomalous coupling to fermions, and the usual MHD equations are fully valid. Any 
source term that was inducing a non-vanishing chemical potential disappears. Thus, the 
transformation of the magnetic fields into fermions is no longer possible. It seems, therefore, 
that the matter fluctuations after the phase transition will be given by the fluctuations right 
before the phase transition. The last statement is in fact wrong for two reasons. First, if the 
phase transition is weakly first order, so that sphaleron processes are in thermal equilibrium 
after it, then any fluctuations of the fermionic charges will disappear. In this particular case 
all anomalous effects that existed before T c are simply "forgotten" , since the system passes 
through an equilibrium period with respect to fermion number non-conservation. Let us 
admit that the electroweak phase transition is strongly first order and a necessary condition 
for EW baryogenesis |36[] is satisfied. Then, there is an important "storage" effect (and we 
come to the second point), which amplifies the estimates of Eqs. (|2.26|) and ( |2.27| ) by many 
orders of magnitude. The point is that the fermion number can sit not only in the fermions 
(and in their associated chemical potential), but also in the hypermagnetic field itself. At the 
EW phase transition, this fermion number must be released in the form of real fermions, just 
because the ordinary magnetic field, which survives after transition, cannot carry fermion 
number. To compute the density of the Chern-Simons number rics of the hypercharge field 
configuration before the EW phase transition we just integrate Ey ■ TCy over the time: 

An cs (t c ) = t E Y ■ Hydt. (2.28) 
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In order to estimate this integral we have to solve the coupled system given by Eqs. (|2.22f) 



and ( |2.24D . The main contribution to this integral comes from the largest time t ~ t c , where 



reactions with right electron chirality flip are in thermal equilibrium. Thus, we can use again 
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the adiabatic approximation (which implies that ~ 0) and obtain: 

AncsiQ = -it: / dt, 

2,71 JO 1 + ljfff c 

A(x,t) = ftyVxfty. (2.29) 

This Chern-Simons number will be released at the EW phase transition in the form of 
fermions, which will not be destroyed by the sphalerons if the phase transition is strongly 
first order. The density of the baryonic number is just given by integrated anomaly: 

n B (t e ) = -^-AncsiQ. (2.30) 

This equation is our main result. Once the hypermagnetic background is specified, the time 
integration appearing in ( [2.29|) can be performed. If the typical scale of the configuration is 
larger than the magnetic diffusivity distance 

r CT ~7^~l(T 9 xL ew , (2.31) 

where L ew ~ 3 cm is the size of the EW horizon at T c ~ 100 GeV, then all the modes of 
the hypermagnetic field with momentum k smaller than 



K ~ Mt, Mo = Mp l ~ 7.1 x 10 17 GeV 
V M o 1.66 jN^f 

will remain frozen into the EW plasma element. Thus the baryon-number fluctuations can 
be written as 

n B \,^, s a' n f H Y ■ V x H Y TM ^ 783 a' 2 \H Y \ 2 



5 — )(x,t c ) = - 1 \ T n = ^r — h^- (2.32) 

V s J 2ncr c s T + T n T£ 22 <r c 7i z 

Notice that in Eq. (pT32|) there is an enhancement by a factor ~ TMq/T 2 arising from the 
time integration of the anomaly term. We also point out that for T w <r the rate of right 
electron chirality flip cancels out. This last expression can be easily written in terms of 
the corresponding curved space quantities and the only point to be kept in mind is that, 
in curved space, the chemical potential is multiplied by the scale factor (which breaks the 
conformal invariance of the AMHD equations). 
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4 Stochastic hypermagnetic backgrounds 



Two qualitatively different classes of hypermagnetic backgrounds can be studied. The first 
class is characterized by a non- vanishing magnetic helicity (i.e. (WyVxWy) 7^ 0), which im- 
plies that the hypercharge field is topologically non-trivial and parity-non-invariant. There- 
fore, in this class of backgrounds not only fluctuations in the baryon number will be produced, 
but also the generation of the baryon asymmetry is possible. We will discuss this possibility 
later in Section 6. 

The aim of this section is to relate the properties of stochastic background with zero 
average magnetic helicity to the baryon number fluctuations. For this type of magnetic 
field {S(ub/ s)(x,t c )) = but {5(riB/s)(x + r, t c )5(ns/s)(x, t c )) 7^ 0, so that only the in- 
homogeneities of baryonic number are produced. We will be interested here in the formal 
aspect of this relation, and we will focus our attention on the study of the phenomenological 
applications in Section 5. 

Consider a stochastic hypermagnetic field whose (parity-invariant) two-point function is, 

Gijif) = (HifflHjix + r)) = J e^Gijifitfk, (4.1) 
where, because of transversality of the magnetic field 

^) ■ (4.2) 

The average appearing in ( |4.1|) denotes an ensemble average. As was previously stated, the 
average hypermagnetic helicity is (A(x, t c )) = (H ■ V x H) = in the case of the transverse 
and parity-invariant two-point function given in Eqs. ( f4.1| ) and ( |4.2| ). 



The assumption of the stochasticity of the background implies that the higher-order 
correlation functions of the magnetic hypercharge fields can be computed in terms of the 
two-point function (|4.1|). For example, the four-point function can be completely expressed 



in terms of the two-point function ([1 . 1|) : 



(Hk(x')Hj (x) Hi (y 1 ) H n (y) ) = [(H^H^x))^)!!^) + 
(H k {x')H l {y'))(H j (x)H n (y)) + (H k (x')H n (y))(H j (x)H l (y')) ]. (4.3) 
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We are now going to compute the level of the induced fluctuations by the above-mentioned 
stochastic hypermagnetic background. We parametrize the spectral properties of our corre- 
lation function by assuming a power law behaviour of its Fourier transform: 



This representation only depends upon two unknown parameters, namely the slope (a) and 
the amplitude, which can be changed by changing k\. The exponential damping appearing 
in the mode function is not the result of any assumption, but it is a direct consequence of the 
fact that, according to the hyperdiffusivity equations ( p. 19] ) and (|2.22|) , all the modes k > k a 
decay thanks to the finite value of the conductivity. The hypermagnetic energy density is 
obtained by tracing the Green function defined in Eq. ( [4.1| ) for r = 0, 



Because of the exponential damping, this integral is always ultraviolet convergent and can 
be very simply performed: 



where £ = k\/T. 

We will often have to compute various four-point functions, and it is sometimes of great 
help to evaluate the higher-order Green functions not in Fourier space (where complicated 
convolutions would appear) but directly in real space. A generic rotationally and parity 
invariant Green function can always be written in real space as 




(4.4) 




(4.5) 




(4.6) 




(4.7) 



where 




d_ 

Or 2 



[r 2 h(r% F 2 (r) 



Or 2 







[h(r% rg{r) 



d r 2 







[r 3 h{r 2 )l 



(4.8) 



and Q{r) is nothing but the trace of our two-point function, namely 




(4.9) 
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(with q — y~ and R — ^-). Clearly, this representation is transverse (i.e. -^jGij(r) = 0), and 
moreover the integral over the spectrum appearing in Eq. ( |4.9| ) can be exactly performed in 
terms of known special functions: 

m = mF (-1-^4) exp [-*] , 5(0) = 4< ( 2 + |) 

(4.10) 

(F(a,b,z) is the confluent hypergeometric function and T(z) is the Euler gamma function 
00)- 

Some physical considerations are now in order. In our problem the relevant scales are 
those are not erased by the plasma conductivity, namely, from Eq. ( |2.31|) , all the scales 
r > r a . Therefore, the physical limit of all our correlation functions will always be the large- 
R limit. Moreover, a physically realistic situation does correspond, in our considerations, to 
the case where the Green functions are decaying at large length scales. If the Green functions 
decay at large distances we automatically exclude the possibility that the energy spectrum 
of the hypermagnetic inhomogeneities will have some peak at large wave-length. The large 
scales (i.e. R > 1) limit of the normalized trace of our Green functions will then be given by 

9(R) = 0lrf + * a R- {a+4) , 9(R) = ||. (4.ii) 

In fc-space the magnetic energy density per logarithmic interval of frequency is defined as 
P~5| , |Tg| as p(k) = dpn/dhik (where pu = ^(\H(x)\ 2 )). Therefore in our case p(k) ~ 
k A (k/ki) a which implies that "blue"( ct> —4) or "violet" (a ^> —4) logarithmic energy 
spectra correspond to the physically interesting case of two-points functions decaying at large 
scales whereas for a < —4 we have "red" logarithmic energy spectra which are connected 
with Green's functions decreasing at small scales. The case of flat logarithmic energy spectra 
(a ~ —4) may quite naturally appear in string cosmological models P~5|. 

It is important to point out that if the Green functions decay at large distances, then 
g(R) < 1. This observation will be of some relevance when we will have to explicitly evaluate 
our fluctuations, since g(R) will turn out to be a useful and natural expansion parameter 
(see Appendix B for details). 

We can also give the large R expression of h{r 2 ) since it can easily be deduced according 
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to Eq. 



p(3\ 2 | e +l 

h(R)=G(0) R-% (4.12) 



where we defined a = —4 + e. 

We are now ready to compute the level of baryon-number fluctuations induced by our 
stochastic background of magnetic hypercharge fields, which is defined by the correlation 
function 

A(r, Te ) = ]j\(s(Pf) &,t e )5(pf) (x + r,t c ))\. (4.13) 

There are two regimes where this calculation can be performed depending upon the relative 
weight of the two rates appearing in Eq. (|2.24 ). If T > the major technical problem we 



have to face is to evaluate the correlation function of the magnetic helicity at two different 
points; this involves, ultimately, the calculation of a four-point function. The algebraic 
details of this long but straightforward calculation are given in Appendix A. The result is 
given in terms of the functions appearing in the real space parametrization of our Green's 
functions given in Eq. 



(H ■ V x H\ (x) (H - V x H) (x + 



A^ / ,dF 2 (r) „ . ,d 2 F 1 (r) „ , ,dF 2 (r) 

-Fi(r) — y—^ ~ 2F X (r) * V ' + 4r 2 [F 2 (r)] 2 + 2rF 1 (r) " 



r dr dr 2 - . ^ r 

-6rF 3 (r)^i + 6Fi(r)F 2 (r) + 2f^i] . (4.14) 
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In this and the following formulae we will often use the notation B{r c ) = H(t c )/T (t ( 
which is convenient since in B the time dependence of the scale factors cancels and the only 
time dependence left is due to the evolution of the effective number of the relativistic degrees 
of freedom in the plasma, N e ff(r). Inserting now Eq. (gjJD into Eq. (gj) we obtain: 

45n f a> T c M ^C(e) h 

A(r ' tc) = ^^^(^ (1 + 0(A)) ' t = T 



C(e) - 



7ie(e + 2) , / T \ 2 , r 



A~ — (-)"-. (4.15) 



T(^) \ (3-e) ' \T n J V, 

Notice that this expression holds for large rate (i.e. F > T-n) and for large scales (i.e. 
g(R)<l,R>l). 
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The second regime in which one may wish to compute the level of induced fluctuations is 
the one where T < E^. The main mathematical problem will be to evaluate the correlation 
function of the hypermagnetic helicity divided by the hypermagnetic energy density. This 
is of course a strongly non-linear object which we cannot compute exactly. By working at 
large scales and in the hypothesis that the Green functions are decaying at large distances, 
the bottom line of this calculation (reported in Appendix B) is : 

B- V x B\ ._. [B-VxB\^ 
— - (x) - [x 

\B\ 2 J V \B\ 2 ) 
((B ■ V x B)(x)(B • V x B)(x + r)} 

(Mm 2 

Notice that to estimate the numerator appearing at the r.h.s. of Eq. ( |4 . 1 , Eq. ( f4.14p can 
be used together with the considerations reported in Appendix A. 



+ 0(g(R)). (4.16) 



5 Phenomenological implications 

Having set all the formalism for computing baryon-number fluctuations, we now come to 
the physical consequences. As we argued in Section 3.2, these fluctuations survive after the 
electroweak phase transition only if it is strongly first-order; we will assume that this is the 
case. We will argue in Section 6.1 that strong enough magnetic fields make an EW phase 
transition strongly first-order even in the case of the minimal standard model. Otherwise, 
some extension of it can be considered. 

An essential quantity entering all expressions for baryon-number fluctuations is the ratio 
between the perturbative and non-perturbative rate of the right electron chirality flip. Fluc- 
tuations are larger for T > T^. In Sections 5.1 and 5.2 we will assume that this is indeed 
the case, and analyse this assumption in detail in Section 5.3. 

As a preliminary warm up let us estimate the amplitude of baryon-number fluctuations 
at the magnetic diffusivity scale for a flat spectrum of magnetic fields (e C 1). If the 
energy sitting in the background magnetic field is comparable with the energy density of the 
photons, (H Y ) ~ ^ci then for the smallest possible scale r ~ l/k a ~ 10~ 9 x(EW horizon 
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~ 3 cm) we get, from Eq. Q4.15 ): 

Air^^^Mr,^. (5.1) 

This estimate is certainly quite large and it is unlikely to be correct, since for such huge 
fluctuation the back-reaction of the created fermions on the magnetic fields and on the dy- 
namics of the electroweak phase transition (which we ignored) must be taken into account. 
Nevertheless, it shows that considerable inhomogeneities in the baryonic number are pos- 
sible on small scales. The estimate ( |5.1| ) considerably exceeds the measure of the baryon 
asymmetry of the Universe n^/s ~ 10~ 10 , thus small size matter-antimatter are possible 
at the EW scale. At the same time, for even larger scales (possibly relevant for structure 
formation), the fluctuations of Eq. ( |4.15| ) are quite minute (since their amplitude decreases 
with the distance as l/r 1+e ) and may be safely neglected. 

The above estimate suggests that a quite natural outcome of the presence of stochastic 
background of the primordial hypercharge field may be a rather inhomogeneous distribution 
of matter and antimatter domains for scales inside the EW horizon. The fluctuations esti- 
mated in Eq. fl4.15| ) are also illustrated in Fig. |1] where the level of fluctuation is plotted 
for different slopes and amplitudes of the hypermagnetic energy spectrum. We clearly see 
that by tuning £ (i.e. by tuning the amplitude of the hypermagnetic energy spectrum) the 
level of induced matter-antimatter fluctuations can be as large as suggested by the estimate 
ofEq. flOp. 

We will now discuss the relevance of the generated fluctuations for the BBN. In fact 
sizable matter-antimatter fluctuations can provide a new type of initial conditions for non- 
homogeneous BBN. From a more conservative point of view, we can instead assume that 
the BBN was essentially homogeneous; then our considerations provide a new bound on 
primordial magnetic fields present at the EW epoch. For completeness we will compare the 
bounds arising from the occurrence of matter-antimatter domains for scales of the order 
of the neutron diffusion distance with the bounds for magnetic fields that are at present 
coherent over much larger scales. 
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Figure 1: We plot the level of fluctuations given in Eq. ( 4.15 ) for different values of ampli- 
tudes (£) and slopes (e) of the hypermagnetic energy spectrum for scales slightly larger than 



the diffusivity scale, which is 10 9 x Lew {Lew is the EW horizon, see also Eq. ( |2.31 )). 



These scales will not be washed out by the finite value of the conductivity, and we can im- 
mediately see that for a sufficiently blue spectrum (e< 1) A ~ 10 3 -10 4 as suggested by the 
estimate of Eq. (|5.1[ ) . Looking at the plots from top to bottom, we see that A decreases for 
increasingly violet spectra (e > 1) and for fixed amplitude (£) of the hypermagnetic energy 
spectrum. We notice that the results illustrated in this plot depend very weakly upon the 
value of the hyperconductivity. In particular, in the present plots, we assumed — 70 as 
fiducial value for the hyperconductivity at the EW scale. 
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5.1 BBN and matter— antimatter fluctuations 



The success of the homogeneous and isotropic nucleosynthesis may impose strong constraints 
upon the baryon-number fluctuations possibly produced prior to the formation of the light 
nuclei. Broadly speaking, the predictions of homogeneous BBN for the primordial abun- 
dances of the light elements are compatible with the observations only if the baryon-to- 
photon ratio lies in a quite narrow range around ns/n^ = 3 x 10~ 10 — 10~ 9 [39 . 

Generally, if A(r,t c ) > — for some length scale r, we have to conclude that matter- 
antimatter domains will be formed. If, on the other hand, A(r, t c ) < — at all scales, only 
positive-definite fluctuations in baryonic density are produced. 

Upper and lower limits on the scales over which a perturbation in baryon number can 
affect nucleosynthesis through neutron-proton segregation are determined by the comoving 
diffusion lengths of neutrons and protons at the beginning of nucleosynthesis. At high tem- 
peratures the diffusion lengths of neutrons and protons are almost the same, since neutron- 
proton equilibrium is guaranteed by weak interactions. After weak interactions have fallen 
out of equilibrium, nucleons retain their identity as neutrons and protons, and diffusive seg- 
regation can occur. Coulomb scattering between protons and electrons (or positrons) give a 
cross section roughly equal to the Thompson cross section. Since neutrons have a magnetic 
moment they scatter electrons with a cross section of 8 x 1CT 31 cm 2 . Neutrons scatter also 
nucleons and the scattering cross section in terms of the triplet and singlet scattering lengths 
is roughly 2.3 x 10 -23 cm 2 . Once the cross-sections of the processes involved are known, the 



diffusion scale is simply given by L{tjj) ~ ^6tuD(tu), where D{t v ) is the diffusion coeffi- 
cient (usually related to the mobility through the Einstein coefficient [^(J) at any given time 
tu = M /T 2 . 

At the onset of nucleosynthesis (T/vs — 0.2m e , where m e is the electron mass) the comov- 
ing diffusion scale turns out to be |^T], g2| 43, 44] 3 x 10 5 cm. The neutron diffusion length 



blue-shifted at the time of the electroweak phase transition is 

L diff (T c ) = 0.3 cm, for T c = 100 GeV . (5.2) 

If A(L^^j,t c ) > — matter-antimatter domains will not be erased by the nucleosynthesis 
time and, at the same time, fluctuations occurring over scales smaller than Ldiff(t c ) at the 
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electroweak epoch are likely to be dissipated ||T], f|2 



Taking again the flat spectrum for magnetic fields and assuming that their energy is 
~ T 4 , we thus obtain for the baryon-number fluctuations at that scale 5(tib/s) ~ 10~ 5 ^> 
10~ 10 . If magnetic fields are large enough (with sufficiently flat spectra), domains of matter 
and antimatter may exist at scales 5 orders of magnitude larger that the neutron diffusion 
length 0. To our best knowledge, there have been no studies of non- homogeneous BBN 
with this type of initial conditions. Of course, there were a lot of investigations of non- 
homogeneous nucleosynthesis, motivated by first-order quark-hadron phase transition [f45|. 
In particular, baryon-number fluctuations with spectral amplitudes growing in frequency 
(and then decaying at large length scales) were recently addressed j43|, with the result 
that these fluctuations are allowed, provided they occur at scales smaller that the neutron 
diffusion length. However, Refs. [41, [El 53, essentially considered positive-definite 



baryon number-fluctuations, rather than with matter-antimatter domains. These results 
were also used in order to constrain the possible baryon-number fluctuations arising in the 



context of topological defects models of baryogenesis ||46|| . It would be very interesting to 
see whether matter-antimatter domains may change BBN bounds on the baryon-to-photon 
ratio by changing the related predictions of the light element abundances. This possible 
analysis will not be attempted here. 

We can instead adopt a more conservative attitude and require that no matter-antimatter 
domains larger than the nucleon diffusion scale exist at the onset of nucleosynthesis. This 
will give some constraints on the primordial hypercharge magnetic fields. In order not to 
have matter-antimatter domains affecting BBN, we therefore demand 

A(L diff ,t c ) < ?j . (5.3) 

From Eq. (|4.15| ), imposing the bound ( |5.3|) , we can translate the constraints coming from 
homogeneity of BBN into an exclusion plot in terms of the only two parameters (e and £) 
characterizing our stochastic hypermagnetic background: 

log£< (-6.262 + l g^ + iloge+(14.88)e + log [Q B h 2 100 ]\ / (A - e) . (5.4) 



3 Note that the energy fluctuations of the electroweak horizon scale are always sufficiently small (i.e. 
Spp/pp <C 1 for r ~ Lew) so that black-holes formation is not expected. 



25 



This condition is reported in Fig. |2| where, for the validity of our approximations, we 
considered the case where e > 0.05 (for e < 0.05 and scales L di ff/r a > 10 9 the corrections 
appearing in Eq. ( |4.15| ) are no under control). We plotted our bounds in the case 0.05 < e < 1. 
There is a second constraint which one might want to impose on our background, namely 
the one coming from the critical energy density: 

p H (t c ) < p 7 (t c ), PH (t c ) = ^(\H(x)\ 2 ), Pl {t c ) = ^AT e// T c 4 (5.5) 

(where (\H(x)\ 2 ) is given by Eq. (|4.5Q . Using now Eq. ( |4.6|) we can convert ( |5.5|) into a 
further (but milder) constraint on our parameter space 

log£ < (^Z 1 -^r + ^°g2 - log^] + loge) /(4- e) (5.6) 

This bound is also reported in Fig. |2] (dashed line). Sizeable matter-antimatter domains are 
produced when the spectrum is sufficiently flat. This feature can also be deduced from the 
level of fluctuations for scales of the order of the neutron diffusion distance, which we plot 
in Fig. |3[ We see that it is quite possible to get A ^> 10~ 10 around the neutron diffusion 
distance. 

We want to notice that an artificial way of relaxing our exclusion plot ( [5.4| ) could be 
(trivially) to enhance the level of the baryon asymmetry by enhancing Qb^qq (up to values 



of the order of 0.1-1). This phenomenon turns out to be similar to the one discussed in [[17 



where it was argued that baryon-number fluctuations with blue frequency spectra might offer 
an interesting mechanism for accounting for large amounts of baryonic dark matter. This 
is purely an analogy, since the problem we are discussing (as we stressed) is not the origin 
of the baryon asymmetry, but the possible bounds on the hypercharge fields. Therefore, 
the value of VtBh\ m is an external parameter for us, but not a computable number (see also 
Section 6). 

A comment is now in order concerning the phenomenological estimates we made in this 
subsection. The system of Eqs. ( 2.22| ) and ( p. 24 ) was solved in the approximation of local 
thermal equilibrium and the possible back-reaction effects were ignored. Numerical solutions 
of this non-linear system of partial differential equations are required if the level of induced 
fluctuations gets too large (i.e. A ~ 1). We will not discuss here how to address this 
complicated numerical problem, but we will come back to it in Section 6. 
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Figure 2: We plot the constraint on the stochastic background of magnetic hypercharge field 
derived in Eq. (|5.4j ) by requiring the homogeneity in the baryon-number fluctuations at the 
neutron diffusion length. We have chosen VtBh\ m = 0.01 and fio^-ioo ~ 1- We notice that by 
changing h wo in the range 0.4 < /i 100 < 1 the quantitative change in the bound is negligible. 
We also took a c /T c ~ 70 as fiducial value. The variation of a c in a plausible range does not 
alter the features of the present plot. The shaded region in the parameter space corresponds 
to matter-antimatter fluctuations that will be erased by the nucleosynthesis time, whereas 
in order to have sizeable matter-antimatter domains at the onset of BBN we have to go 
above the line. In the dashed line we report the critical energy density bound given in Eqs. 



D and (O) 
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Figure 3: We plot the expected level of fluctuations given in Eq. Q4.15Q for scales of the order 
of (and larger than) the neutron diffusion distance given in Eq. ( |5.2j ). In the plot r n = L^//. 
In this case we always took a c /T c ~ 70. We can notice that for flat enough spectra a very 
interesting level of fluctuations is allowed. For violet spectra the fluctuations are certainly 
suppressed. The same trend exists for even larger scales. 
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5.2 Comparison with other bounds on primordial magnetic fields 

In this section we are going to compare the BBN bound reported in Eq. ( |5.4|) with other 
possible bounds, which could be applied on primordial magnetic fields. 

It is well known that there are direct bounds on primordial magnetic fields at the nucle- 
osynthesis epoch [p0[]. Moreover, quite recently, two constraints on magnetic field intensities 



were derived using, respectively, the anisotropies in the microwave sky |2T[ and the Faraday 
rotation correlations ||24j| . These last two bounds apply of course to the case of fields that are, 
today, coherent over length scales much broader than the (present) nucleosynthesis scale. In 
order to see how stringent our bound is, we should compare it with the ones already available 
and coming from larger scales. In this sense our aim is to show that our bounds are clearly 
more stringent for small-scale fields but cannot compete (at even larger scales) with the ones 
coming from the CMBR anisotropies and from the Faraday rotation measurements. 

Since the CMBR is isotropic to a very high degree of accuracy, its small anisotropies can 
constrain the intensity of a constant magnetic field (coherent over the present horizon size 
2lfl), which could modify the evolution equations of the matter sources by introducing a 



slightly anisotropic pressure fl22|. The calculation of the CMBR anisotropies can be carried 
out also in the case of slightly skew stresses, whose numerical weight depends upon the 
magnetic field intensity. By comparing the final result with the level of anisotropies detected 
by COBE, it is possible to compute how big the magnetic field intensity should be not to 
conflict with observed anisotropies in the microwave sky. At the present time, the constraint 
is \H (t )\ < 6.8 x 10~ 9 (f2 ^ioo) 1//2 g auss over a length scale L (t ) ~ 9.25 x 10 27 ft^o cm - 
The authors of Refs. |21], |22| gave the bound in terms of h 5Q (the present uncertainty 
in the Hubble parameter in units of 50 km/Mpc sec). For consistency with our notation 
and in order to make the comparison with other bounds easier, use instead h wo = h 50 /2 
taking, as usual, 0.4 < h 1Q0 < 1. By blue-shifting the bound of up to T c we get 
\H {t c )\ < 1.12 x 10 22 (fio^ioo)^ 2 gauss at a scale L (t c ) ~ 7.19 x 10 12 /ij~ o cm > where we 
assumed that the magnetic field scales as l/a 2 (r), as it is plausible to demand for length 
scales larger than the magnetic diffusivity scale. We also took into account the evolution in 
the effective number of degrees of freedom in the plasma, which is N e ff(t c ) ~ 106.75 and 



N e ff(tdec) ~ 3.90. The bound on H (t c ) turns into a bound on B(t c ) = J(H 2 )/T 2 (recall 
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Figure 4: We plot our bound (lower curve) obtained in Eq. (|5.4|) and the bound derived in 



(upper curve) from the CMBR anisotropies in the case of a magnetic field coherent 
over the horizon size today. For values of e and £ above the shaded region the level of 
fluctuations exceeds the bound ( |5.3| ). Again the variation of VL B h\ m (possibly between 0.1 
and 0.01) changes the plot by only a few percent. In this plot we took h\w = 0.6. 

that 1 gauss = 1.95 x 10- 20 GeV 2 ). 

From Eq. ( |4.6j ) (taking r ~ L (t c ), T ~ T c ) it actually results, that the parameter space 
of our model has to satisfy the following requirement 

log£ < (-2.30 + l - log [n h 2 100 ] + X - loge + [l4.13 - 1 

For comparison this constraint is reported in Fig.|I], together with the bound of Fig.|2|. Since 
the region defined by Eq. ( |5.4| ) is always below the curve of Eq. (|5.7|), we conclude that 
the bound imposed by the homogeneous BBN is more constraining than the one reported 
in 



log h 



100 



(5.7) 



21, B2fl for e > 0.05. Clearly, for some very small slopes, the COBE bound will become 



better, but we cannot compute this critical value of e since our approximation breaks down 
before it. In P3 the polarized emission of few hundred galaxies was reviewed and another 
bound on the present intensity of large-scale magnetic fields was reported. A large-scale 
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field should produce an additional shift in the polarization plane of the incoming radiation 
according to the Faraday effect. In fact the polarization plane of the synchrotron radiation 
gets shifted, in the background of a magnetic field, by an amount that is directly proportional 
to the integral of the magnetic field times the electron density along the line of sight. By 
subtracting, from the total angular shift, the one produced by each galactic field, it is possible 
to constrain the intensity of any field coherent over scales larger than the galactic one. The 
only uncertainty with this procedure is that the measurements assume that the magnetic 
fields of the Milky Way and of the other galaxies are known to a very high degree of accuracy, 
since they have to be subtracted from each estimate of the Faraday rotation. The constraint 
obtained with this method turns out to be {H^to)] < 1 X [ 2 - 6xl0 _ ?cm 3 }h W0 x 10" 9 gauss, 



II B 

-1 



for fields now coherent over scales Li(to) — (10 — 50)/ij~ o Mpc. Assuming that the mean 
(present) baryon density is n B ~ 1.13 x 10~ 5 (Q B h1 00 ) cm" 3 , the (blue-shifted) field intensity 
and its coherence length will be respectively \Hi(t c )\ < 3.79 x 10 19 (fi,B/iioo) _1 ^ioo gauss 
and L x {t c ) ~6x 10 10 ^roo 

cm where we took ( Lo(^o) ~ ^ Mpc). By translating this bound 
in the (£, e) plane, we obtain the following relation 

log£< (-4.77 + l - loge -^logQ B /4 +[l3.09-^ log h m ] e ) / ( 2 ~ |) • ( 5 - 8 ) 

This constraint is illustrated in Fig. [5], where we see that our curve is always below the 
curve reported in Eq. ( |5.8|) . Our constraint is again more stringent than the one of Eq. 
Q5.8D for e > 0.05. The bounds on ordinary magnetic fields at the nucleosynthesis epoch 
also apply in our case. In order not to affect the universe expansion at nucleosynthesis it 
should hold (see for instance Kernan et al. in Ref. [pOfl ) that pn < 0.27 p v (where pu is 
the magnetic energy density defined in Eq. Q5.5Q ) and p v is the energy density contributed 
by the standard three light neutrinos for T C 1 MeV). Therefore in terms of £ and e this 
bound reads: 

(11.30-ilog%)e + loge-0.2 

log£ < —s . 5.9 

4 — e 

This bound is reported in Fig ^ and is compared with our bound. We see that the bound 
coming from the absence of matter-antimatter domains at the nucleosynthesis epoch is more 
constraining for (by two orders of magnitude for logarithmic energy spectra with e<l). 
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Figure 5: We compare the bounds of Eqs. ( |5.4j ) and ( |5.8| ). The most constraining bound 
is the one given by the lower curve given by Eq. ( |5.4j ). The upper curve derived in the 
context of Faraday rotation measurements is more constraining than the one obtained from 
the isotropy of the CMBR. The numerical value of the parameters for which this plot is 
obtained is the same as for Fig. [2] and [|. 



32 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

e 



Figure 6: We compare the direct bound on primordial magnetic fields at the nucleosynthesis 
epoch derived in Eq. ( |5.9|) with the bound derived in Eq. Q5.4D , which applies to primordial 
magnetic fields. We see that our requirement is again more constraining than the one given 



in Refs. |20] for blue spectra, whereas for e > 1.4 the bound given by ( |5.9| ) is better 
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Figure 7: We plot the constraint derived on the basis of BBN considerations (Eq. ( |5.4|) ) 
together with the requirement derived in Eq. ( |5.10D . The dashed line corresponds to a 
magnetic field of EW origin strong enough to rotate the polarization plane of the CMBR 
under the assumption that the CMBR is polarized. We can see that if such a field has 
an electroweak origin, then also sizeable matter-antimatter fluctuations will be present for 
scales of the order of the neutron diffusion distance. 



Note also that, according to |j48[ , the bound (|5.9| ) may in fact be absent, because there are 
other mechanisms, besides magnetic diffusivity, that can dilute the magnetic fields before the 
BBN. Then our bound remains the only one that can be applied to the small-scale magnetic 
fields. 

Another interesting numerical value of the magnetic field, which might be compared 
with our considerations is\H 3 (t dec )\ ~ 10~ 3 gauss coherent over a scale L 3 (t dec ) ~ 1.690 x 
10 23 (f2o^ioo) -1/ ' 2 cm ' which is the size of the horizon at the decoupling. If this field would 
be present at the decoupling time (when the temperature was roughly Td ec = 0.258 ev) 
it might also rotate the polarization plane of the CMBR provided the CMBR is weakly 
polarized [p3fl . At the scale of the electroweak phase transition the blue-shifted field and 
the corresponding correlation length read respectively \H 3 (t c )\ ~ 1.36 x 10 21 gauss and 
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L 3 (t c ) ~ 1.447 x 10 11 (Slo^ioo) ^ 2 cm ano - this imposes for our parameters the following 
requirement 

log£> (-3.22 + i log e + 13.28 e - | log Moo]) / ( 2 " ■ ( 5 - 10 ) 

It is of some interest to notice from Fig. [J] that the region defined by Eq. ( p.lCj ) falls 
in the forbidden area of Fig. ^. This means that the BBN bound of Eq. ( |5.4j ) excludes 
the possibility that a primordial magnetic field of EW origin is strong enough to rotate the 
polarization plane of the CMBR. On the other hand, if non-standard initial conditions for 
the inhomogeneous BBN scenario (i.e. matter-antimatter domains) would be allowed, this 
conclusion might be relaxed and the existence of such an intense field at the decoupling 
epoch might be accommodated. It is anyway amusing that in our present discussion the 
existence of a magnetic field at the decoupling epoch might imply the presence of small-scale 
antimatter domains at the onset of BBN. 



5.3 The rate of right-electron chirality flip 

In Section 3 we pointed out how important is, in our context, the interplay between the "per- 
turbative" rate given by the right electron chirality flip processes and the "non-perturbative" 
one coming from the anomaly. For the MSM the perturbative rate of chirality flip has been 
computed in and is determined by the right electron Yukawa coupling. If MSM is a 
correct theory, then T > Y-h only for extremely small magnetic fields, H 2 /T 4: < j§^^j^ — 
10~ n . So weak hypermagnetic fields do not produce any interesting fluctuations. For larger 
magnetic fields the approximation outlined in Eq. Q2.27 ) must be used. The amplification 
factor that appears in Eq. (|2.32|) can be extracted from f3l| : 



f tc , TM /100 GeV\ r/ . 

J dtT~^-~ 350 [ J [(-1.1 + 2.4 z H ) + 1.0 + fc t 2 (0.6- 0.09 (5.11) 

where xh = limr-*oo mf ^ T ' > is the high-temperature limit of the Higgs thermal mass and h t 
is the top-quark Yukawa coupling. Taking xu ~ 0.6 (for which the scattering contribution 
to the rate is always dominant with resepct to the decay contribution), we find that the 
integrated rate is 655 x (100 GeV/T c ). With the use of this number the analysis of Section 
5 can be redone with the result that no interesting baryon-number fluctuations can be 
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produced from stochastic hypermagnetic background. So, for MSM, one hardly expects any 
cosmological consequences coming from the background of the type ( f4.4p (for other types of 
primordial hypermagnetic fields, considered below in Section 6.2, the conclusion is different). 

However, in the extensions of the standard model, the rate T can be naturally larger than 
in the MSM. For example, in the context of the Minimal Supersymmetric Standard Model 
(MSSM) the right-electrons Yukawa coupling is enhanced by a factor 1/ cos j3 (tan/5 gives 
the ratio of the expectation values of the two doublets), so that Tr can be larger by a factor 
10 3 for experimentally allowed values of tan/3 ~ 50. Moreover the right-electron number is 
now shared between electrons and selectrons, and it is necessary to consider also processes 
that change the selectron number. 

The question we now want to address is more phenomenological. Namely we want to 
see how large the perturbative chirality flip rate T should be in order to produce sizeable 
matter-antimatter fluctuations, which could influence the BBN. For this purpose, we just 
require that T > T^, with taken from Eqs. (|2.32j ) and (|4.6j ), and use the minimal 
amplitude of the hypermagnetic field obtainable from the bound ( |5.4|) , which can produce 
sizeable matter-antimatter fluctuations. This gives T/T c > 10~ 9 , which corresponds to the 
right electron perturbative freezing temperature Tr = 10 5 TeV. As we discussed above, these 
values of the temperature is perfectly possible, say, in the MSSM. If the actual value of the 
freezing temperature is smaller than 10 5 TeV, the stochastic hypercharge background of type 
(O) produces baryon- number fluctuations too small to affect BBN. It is then interesting 
that a quite energetic stochastic hypermagnetic background can be accommodated in the 
MSM without any (potentially dangerous) implications. The energy density stored in this 
background can then be able to influence the dynamics of the EW phase transition without 
conflicting with any bound derived from BBN. This will be one of the subjects of the following 
section. 



6 EW phase transition and baryogenesis 

The aim of this section is the discussion of the influence of the hypercharge magnetic field 
on the electroweak baryogenesis (for reviews, see |^9[). First, we will consider the EW phase 
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transition in the presence of the hypermagnetic field. Then, we will show that the occurrence 
of some specific hypermagnetic configurations in the EW plasma could be responsible for 
the baryon asymmetry of the Universe (BAU). 



6.1 EW phase transition 

The hypercharge magnetic field, present at high temperature, can influence the dynamics of 
the phase transition. The physical picture of this phenomenon is exactly the same as the 
macroscopic description of superconductors in the presence of an external magnetic field. The 
normal-superconducting phase transition, being of second order in the absence of magnetic 
fields, becomes of first order if a magnetic field is present. The reason for this is the Meissner 
effect, i.e. the fact that the magnetic field cannot propagate inside a superconducting cavity, 
and, therefore, creates an extra pressure acting on the normal- superconducting boundary 
50fl . Our consideration below explores this simple picture. 

Consider the plain domain wall that separates broken and symmetric phase at some 
temperature T, in the presence of a uniform hypercharge magnetic field y^. Far from the 
domain wall, in the symmetric phase, the non-Abelian SU(2) field strength (W|) is equal 
to zero, because of a non-perturbative mass gap generation. Inside the broken phase, the 
massive Zj combination of 3^' and W|, 

Zj = cos 6 w Wf - sin Q w y i (6.1) 

must be equal to zero, while the massless combination, corresponding to photon Aj m , sur- 
vives. The matching of the fields on the boundary gives Aj m = yj cos 9w- Thus, an extra 
pressure \\Hy | 2 sin 2 9w acts on the domain wall from the symmetric side. At the critical 
temperature it must be compensated by the vacuum pressure of the scalar field. If we ne- 
glect loop corrections associated with the presence of magnetic fields, then the condition 
that determines the critical temperature is: 

~\H Y \ 2 sm 2 9 w = V(0,T c )-V(<p min ,T c ) , (6.2) 

where V((p, T) is the effective potential in the absence of magnetic field, (p m in is the location 
of the minimum of the potential at temperature T. 
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Figure 8: We plot the requirement obtained in Eq. (|6.2|) in the case of a stochastic hy- 
permagnetic background. The dashed line refers to the case of tuh = 80 GeV, whereas 
the dot-dashed line refers to the case of m# = 160 GeV. With the full line is reported for 
comparison the bound coming from BBN and discussed in Eq. ( |5.4| ). We can clearly see 
(shaded region) that for steep enough hypermagnetic energy spectra (i.e. e > 0.4 -0.6) it is 
possible to have a strongly first-order EW phase transition consistent with the bound (|5.4|) . 

The above consideration was dealing with the uniform magnetic fields. Clearly, it remains 
valid when the typical distance scale of inhomogeneities of magnetic field are larger than the 
typical bubble size. This is the case for bubbles smaller than the magnetic diffusion scale, 
and, in particular, at the onset of the bubble nucleation. Thus, the estimate of the critical 
temperature coming from ( |6.2|) is applicable. For bubbles larger than the diffusivity scale, 
the presence of a stochastic magnetic field will considerably modify their evolution. In 
particular, the spherical form of the bubbles is very likely to be spoiled. 

Relation ( |6.2| ) may be used to define a minimum magnetic field, which can make a 
phase transition strong enough to allow electroweak baryogenesis in the MSM. One can just 
fix the Higgs mass, find the temperature at which the minimum of the effective potential 
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satisfies the constraint (p min /T > 1 |36]], and read off the hypermagnetic field from Eq. (|6.2| ). 
With the use of the two-loop effective potential computed in |5"T| , we get, for m H ~ 80 GeV, 
(\H Y \ 2 )/T* > 0.06, whereas if m H ~ 160 GeV we will have (\H Y \ 2 )/T* > 0.3. For stochastic 
backgrounds these constraints are plotted in Fig. |8] in terms of our variables £ and e, 
characterizing the spectrum. We see that the region of parameter space where a strongly 
first-order EW phase transition is possible, without spoiling BBN with excessive matter- 
antimatter domains extends, from e ~ 0.5 up to e > 1. Therefore we come to the conclusion 
that for violet hypermagnetic energy spectra the level of induced fluctuations is quite tiny 
at the neutron diffusion distance, but the dynamics of the phase transition can be strongly 
affected. The magnetic fields, which can modify the nature of the phase transitions do not 
appear to be subjected to any other constraints. 

The observation that the presence of primordial hypercharge magnetic fields at the elec- 
troweak epoch may make an EW phase transition strongly first-order removes a main ob- 



jection against the possibility of baryon asymmetry generation in the MSM [3g, (49[ . We 
also note, that the background magnetic field breaks C and CP symmetries, which may 
considerably change the analysis of different processes near the domain wall, which are used 
in EW baryogenesis mechanisms. This study is beyond the scope of this paper. Instead, we 
will point out in the next subsection that some specific configurations of the hypercharge 
magnetic field may create the net baryon asymmetry. The following discussion is similar to 
considerations of the baryogenesis from Chern-Simons condensate in |3"6] . 



6.2 Chern-Simons condensates and the BAU 

In Sections 4 and 5 we were concerned with the case of stochastically distributed hypermag- 
netic fields. In this case, the average baryon asymmetry vanishes (i.e. (5(nB/s)(x,t c )) = 0) 
in spite of the fact that the fluctuations of the baryon number may be considerable. Thus, 
we assumed that the source of the baryon asymmetry has no relation to the primordial 
hypermagnetic fields. For example, the BAU may have been generated because of GUT 
interactions or at the EW phase transition. 

In this subsection we are going to discuss hypermagnetic backgrounds that may give rise 
to the BAU. We have no idea whether these types of background can or cannot be generated 
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by some mechanism. Our aim is to point out the essential features of the hypermagnetic 
field that are necessary for the production of net baryon asymmetry. 

If the hypermagnetic field configuration is topologically non-trivial (i.e. Hy-Vx^y ^ 0), 
then S(ub/s) 7^ 0. As a particular example we will discuss the case of a Chern-Simons wave 
configuration 

34 = y(t) sin k z, 
yy = y(t) cos k z, 

y z = 0. (6.3) 
The hypermagnetic field is 7i = V x y and the magnetic helicity is simply 

H-VxH,= k H 2 (t) 



with lH.it) = k y(t). Thus from (|6.3| ) we obtain that: 



Hl + Hl + ^ = W(t), n cs = -£- 2 H ■ V = ^-« 2 W 

We notice that this configuration describes a magnetic knot with uniform energy and Chern- 
Simons density over the whole observable Universe. Similar configurations are used in the 
MHD treatment of the dynamo instability 0. Inserting the configuration ( |6.3|) into the 
evolution equations (|2.22| ) and into ( |2.24| ) , we get a system of non-linear ordinary differential 
equations 



(where w = tT). Equations (|6.4|) are in fact known as a generalized Lotka-Volterra sys- 
tem ||52|| . They can be solved numerically for different types of initial conditions. As in 
the discussion of the stochastic backgrounds, we will consider our system in the adiabatic 
approximation and we will then use the general set of equations (|2.22j) and (p.24[) , valid for 
arbitrary backgrounds. 
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As usual, the magnetic diffusivity a c defines the diffusion scale and therefore a Chern- 
Simons wave configuration with typical momentum (fc ) larger than k a ~ 10~ 7 T C will be 
washed out. For smaller ho, we get from ( |2.32| ): 

-Wy s ^m^aW( lt )„»»W«, for r>r K 



s J 2 7r a c \ s ) \T C J \T C J ' \T C J \ T„ 4 

n s\ ± \ . . Utt fn f \ 7 _ n n {T R \ { k 



H 783^lfJ fc » rM ^°-H7fJ^J- for (6 ' 5) 

Let us now assume that we work only in the framework of the MSM. Then, in order to produce 
baryon asymmetry we need a strong enough first-order phase transition, and therefore, a 
strong enough magnetic field (see previous subsection). Thus, T-h > T. Using the fact that, 
in the MSM, T R ~ 80 TeV we see that fcf \ ~ 1(T 10 for k /T c ~ 1(T 12 . This value is well 
below the magnetic diffusivity scale and, therefore, this type of configuration, if ever created, 
will survive till the EW epoch. 

In the extensions of the standard model one may have a strong enough electroweak phase 



transition without any magnetic field [|53]]. In addition, the perturbative electron chirality 
rate may be considerably higher than in the MSM. Thus, the hierarchy < T may be 
naturally realized. Then, for f |?-J \ — 10~ 20 the BAU calculated from the hypermagnetic 
field is of the right order of magnitude. For example if the typical momentum fc of the Chern- 
Simons condensate is of the order of the EW horizon at T ~ T c (i.e. ko ~ ~ 10 _16 T C ) 
then for small enough hypermagnetic energy (i.e. 7i 2 — 10~ 4 T C 4 ) the BAU is ~ 10~ 10 . 
Thus, it is not excluded that the baryon asymmetry of the Universe is a consequence of the 
topologically non-trivial primordial hypercharge magnetic field. 



7 Concluding remarks 

There are no compelling theoretical reasons against the existence of long lived Abelian hy- 
percharge fields at the electroweak epoch. In the present paper we showed that, if they did 
exist, they could have a number of cosmological consequences. The stochastic hypermagnetic 
backgrounds induce baryon-number fluctuations because of the electroweak anomaly. These 
fluctuations may produce sizable matter-antimatter domains at the onset of BBN and affect 
its predictions. Magnetic fields can change considerably the dynamics of the electroweak 
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phase transition in the MSM and make it strongly first-order even for large Higgs masses. 
Topologically non-trivial hypermagnetic configurations may be responsible for the BAU. 

We left aside a number of questions that may be subjects of future studies. For example, 
in the study of AMHD equations, we focused our attention on the case where the correlation 
scale of the velocity field was much smaller than that of the magnetic field, and we also 
assumed that the velocity field was (globally) invariant under parity transformations (i.e. in 
the absence of global vorticity over the EW horizon, at the epoch of the phase transition). 
Owing to different phenomena (e.g. bubble collision) turbulence may arise inside the EW 
horizon, leading to a non-mirror symmetric velocity field over some length scale typical of the 
mechanism responsible for the turbulence. If the turbulence produces a non-zero vorticity of 
the plasma, then the collective plasma motions may be transformed into fermion number via 
the amplified hypermagnetic field through a kind of EW dynamo mechanism. We completely 
neglected the possible occurrence of (global) vorticity, and to relax this hypothesis may be 
interesting. 

We do not know what is the possible influence of matter-antimatter domains on the 
inhomogeneous BBN scenario. In particular we have no idea if some spectral distribution 
of hypermagnetic fields could induce a spectrum of baryon-number fluctuations, which can 
lead to a larger baryon-to-photon ratio. 
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A Magnetic helicity correlations (T > T#) 

The aim of this Appendix is to compute explicitly the two-point correlation function of the 
magnetic helicity A(x), namely 

(A(£)A(£ + r)) = ((H ■ V x H)(x)(H ■ V x H)(x + f)) (A.l) 
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in terms of the two-point function 

G ij (r) = F 1 (r)5 ij + r i r j F 2 (r). (A.2) 

The results illustrated here are quite relevant for a precise estimate of the level of fluctuations 
induced by a stochastic background of hypermagnetic fields. 

The estimate of the correlation function (|A.1| ) may be carried out either in real space 
or in Fourier space. In Fourier space the calculation can be reduced to the estimate of 
a convolution, whereas in real space the main algebraic task is to compute the various 
derivatives appearing in the ensemble average. 

The stochastic average appearing in Eq. (|A.l ) can be rewritten as 



((H ■ V x H)(x)(H ■ V x H){y)) = 
d d 

eijk^mni lim lim {H k (x)H j (x)Hi(y)H n (y)). (A.3) 

x'^xy'^y OX Oy m 

This expression turns out to be quite useful, since it allows us to perform the derivations 
with respect to x % and y m after the average is taken. 

If we now use the fact that the background of hypercharge fields is stochastic, we can 
write that 

{H k (x')Hj (x) Hi (y)H n (y) ) = [ (H k (x')Hj (x))(Hi(y)H n (y) ) + 
(Hk(x*)Hi(tf))(Hj(x)H n (y)) + (H k (x*)H n (y)){H j (x)H l (y > ))}. (A.4) 



Inserting in Eq. (|A.4 ) the representation ( |A.2| ), we get the following relation 



where 



((H ■ V X H){x){H ■ V x H){y)) = 
d d 

eijkemni |im^un— — x ; y, y) + B(x , x ; y, y) + C(x , x ; y, y)\, (A.5) 



A(x,x;y,y) = [F 2 (\& - x \)F 2 {\y l - y\)(x' k - x k )(x'j - xj)^ - y t )(y' n - y n ) 

+ f^x 1 - x\)f 2 {W - y\)5 kJ (y'i - yiWn - yn) 

+ F 2 (|f - xDFidy - y\)5 n i(x' k - x k )(rfj - xj) 
+ F 1 (\x'-x\)F 1 (\y'-y\)6 kj 5 ln }. 
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B(x,x'; y, f) = \F 2 (\x' - jf\)F 2 (\x-y\)(x' k - y' k ){x\ - yfi(xj - yj)(x n - y n ) 

+ Fxdtf - y > \)F 2 (\x - y\)Su(xj - yj)(x n - y n ) 

+ F 2 (|f - #\)Fi{\2 - y\)5 nj (x' k - y' k ){x\ - y[) 

+ F 1 (\x'-y'\)F 1 (\x-y\)5 kl 5 jn }, 

C(x,x';y,g) = [F 2 (|f - y\)F 2 (\x - y\)(x' n - y' n )(x' k - y k )( Xj - y' 3 ){x x - y[) 

+ F 1 {\x > - y\)F 2 (\x-y'\)6 kn (x j - y' 3 ){x x - y\) 

+ F 2 (\x - y'\)F 1 (\x' - y\)5 t j« ~ Vk)K ~ V*) 

+ F 1 (\x*-y\)F 1 (\x-y'\)8 k J jl }. (A.6) 

Recall that, for a generic function f(r) (where r = \ f\, r a = x a — y a ), the following trivial 
relation holds 

d d „. , „ ldf(r) r i r m Qf( r \ r l r m d 2 f(r) 



After having performed the derivatives in Eq. ( |A.5| ) we can contract the various Levi-Civita 
tensors with the Kroeneker symbols. By then taking the limits indicated in ( |A.3| ) and ( |A.5| ) 
we obtain 

' 'H - V x H) (x) (H - V x H) (x + r 



4^ /N dF 2 (r) „ . ,d 2 F 1 (r) A 2rr , , , l2 „.,dF 2 (r) 

-6rF 2 (r)^^ + 6Fi(r)F 2 (r) + *[^&\ ■ (A.7) 
ar \ dr ) 

In this form the four-point correlation function is completely expressed in terms of the two- 
point function. Of course we stress that this decomposition holds provided the fields are 
stochastically distributed, namely if and only if ( |A.4|) is satisfied. 



B Magnetic helicity correlations (T < T#) 

If r < r# the correlation function appearing in the final expression of the level of the 
fluctuations turns out to be the stochastic average of a quantity that contains the magnetic 
helicity in the numerator and the hypermagnetic energy density in the denominator. Even 
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if this case turns out to be less relevant for the phenomenological considerations presented 
in the main discussion, we want to outline the main idea that can be used in order to get a 
large-scale estimate of 

-jp — )M[ — jp — iMh (B.l, 



The strategy we will use will be to express ( |B.1| ) in terms of an appropriate path integral 



whose functional derivatives will reproduce the correlation function we want to compute. 
From Eq. ( |B . 1|) we have, formally 



H • V x H\ , . i H-V x H 



x l TTo \?2, 



H 2 I \ H 2 

d d ( $ & & $ \ 
lim lim lim lim 6n k 6 abc —^- — \ r ^,^, . r ^,^. r ^,^,. r ^,^. W[ J] (B.2) 
S' 1 -,x 1 s' 2 -*x 2 a^o/3-,o tJk abc dx\ dx a 2 \5J k {$\) 8Jj(xi) 5 J c (x' 2 ) SJ b (x 2 ) J j =0 



where 



r d 3 p 


1 d 3 q J 







W[J\ 

exp{ - l -J d 3 x J d 3 yH i {x)[K{x,y)] lj H J (y) + i J Ji(x)Hi(x)d 3 x 
+ ipi J F;(f)5 (3) (£ - x l )d 3 x + iqj J Hj(x)5^\x - x' 2 )d 3 x 
- a\p\ 2 -P\q] 2 }- (B.3) 

In Eq. ( p.3|) we used the fact that formally holds the following relation 



a^o J r> 



\H\ 2 Atiol^oJ p 

By appropriately redefining the source term in the path integral, Eq. (|B.3| ) can also be 
written as : 



exp{- l -Jd 3 x J d'yH^K^^H^y) - a\pf - f3\qf 

+i J S i (x)H i (x)^d 3 x (B.4) 
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where, in the present case, 

Si(x) = Ji(x) +Pt 5^(x~x' 1 )+q t 5^(x-x' 2 ). (B.5) 

By defining 

hi(x) = J d 3 y[K^{x,y)} i3 H j {y) 

we obtain that the argument of the exponential can be expressed as 

% - J d 3 x{h t (x) - J d 3 xS l {y)[K^{x,y)] l ^ ~\j d " x j d'yS^G^x-yDSM. (B.6) 

Notice that the symbol [K(x,y)]ij must have the following properties, which will be 
important also for the calculation of the functional integral: 

/ d 3 y[K^(x, y)]ij[K*(y, z)\ jk = [K(x, z)} ik 
J d 3 y[K^ (x, y)] tJ [K-± (y, z)] jk = 6*6® (x - z) 
J cPylK-i&yflijlK-ifaz)]^ = [K(x, z)~ l ] lk 

[K- 1 fr$] ij = -G ij (\x-y\). (B.7) 

We can integrate the part that is quadratic in the fields; then W[J] becomes 

d 3 p f d 3 q 



W[J] = W[0]J^J^ex P {-yd 3 xJd 3 yS t (x)G t] (\x-y\)SM 

- «|Pl 2 -/3|g| 2 }, " (B.8) 



where W[0] is the usual Jacobian. Using the definition of Si(x) we obtain for W[J]: 

W\A = f £p f f0_ C( qi , Pi ; qj , Pj ) 

W[0] J p J q 

x exp j-1 J d 3 x J cPyJifflGijQZ - y\)Jj(y) - i J d 3 xJ l {x)L l {x)^ (B.9) 

with 

Li(x) = p m G m i(\x - x[\) + q m G m i(\x - 41) (B.10) 

and 

C{qi,Pi;qj,Pj) = ~qiGij(0)qj ~ ^PiGij{0)pj - qiG^x^ - & 2 \)Pj ~ a \P\ 2 ~ P\q\ 2 ( B - n ) 
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(Grjj(O) comes because there are two delta functions centred at the same point for the terms 
quadratic in q^ and p^). Performing the functional derivatives we obtain (evaluating the 
generating function for J = 0) 

\6J k (x" 1 )6J j (x 1 )6J c (x' 2 )6J b (x 2 ) L Vj =0 
{[Gkj(\xi - xi \)G cb (\x 2 -x 2 \) 

+ G kc (\x[ - x' 2 \)G jb (\xi - x 2 \) + GkbQx'i - x 2 \)G jc (\xi - 

- L h {x 2 )L c {^G ih (\S 1 -x[\) - L fe (f' 1 )L j (fi)G'6 C (|f2 - x' 2 \) 

- L^x^Lk^Gbjdx! - x 2 \) - L h (xi)L k (tf{)G C j(\x-i - x 2 \) 

- L c (x , 2 )L j (xi)G b k(\x'i - x 2 \) - L b (x 2 )L k (tf l )G k c(\x'i - x 2 \) 

+ L fe (x , 1 )L,(f 1 )L 6 (4)L c (x 2 )}W[0]. (B.12) 



Notice that the fifth and sixth terms of Eq. ( |B.12|) vanish when contracted with the epsilon 



tensors. In order to perform the integration over p and q, we have to expand the expression 
giving us the correlation function for 

9(r) = g < 1, (B.13) 

(see also Eq. Q4.ll]) ). This approximation holds for sufficiently large scales, provided the 
Green functions decay for R 3> 1. This requirement is automatically satisfied in our case, 
since we only consider the situation where the energy spectrum is increasing in frequency 
(i.e. blue or violet spectra). We now take the limits for x[ — > X\, and for x' 2 — > x 2 and Eq. 
H2[) becomes: 

( x ) Jp I ( x + r ) 



lim lim-^ fLl^ir) + q' 2 A 2 (r) - p 2 A :) (r) 



d 3 p r d 3 q , 
»3^o 167T 2 J p J q 

+ 2A 4 (r)(g-^]exp{-i^^(g 2 +p 2 )-a|gf-/3|pl 2 + 0(^(r))}, (B.I I.) 



where 

Ai(r) = /((#• V x H){x){H- V x #)(£ + 
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A 2 (r) = A 3 (r)= 1 -[(H\x))} 3 9(r) 

Mr) = ^[(H\x))] 3 9(r). (B.15) 

We can now integrate over q and p. We notice that since the integral in convergent also 
for a, (3 — » the limits can be taken before the integration. It is convenient to perform the 
integration over q and p separately; in this way, after angular integration, the apparently 
Gaussian integrals can be expressed as ordinary exponential integrals of the type 

r A „ e -^ = r( n + i) 

Jo a n+1 v ' 

After integration, A 2 (r) and A 3 (r) cancel whereas the contribution of the term containing 
A 4 (r) vanishes because of the angular integration. The final result obtained in the assump- 
tion that the Green functions decay at large distances is then 

H* ) {X) [ IP J ^ + r )/- 

which is exactly what we report in Section 4. The method used in the present Appendix can 
also be exploited in order to compute further corrections, if needed. 
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